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Abstract 

We describe an extension to the Geant4 software package that aUows it to be 
used as a general purpose X-ray tracing package. We demonstrate the use of 
our extension by building a model of the X-ray optics of the X-ray observatory 
XMM-Newton, calculating its effective area, and comparing the results with the 
published calibration curves. 
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1. Introduction 

The Geant4 simulation toolkit is a software library for high-energy physics 
developed by an international collaboration under the lead of CERN [l| . Aiming 
at supporting the simulation and analysis of data collected by the Large Hadron 
Collider, Geant4 has found a keen following in many other areas where the ability 
to model complex geometries and the interaction of particles and radiation with 
matter plays a prominent role. The open nature of the toolkit, and the large 
libraries of geometrical models and physical processes that are available make 
it easy to use and extend. 

The importance of detailed ray-tracing simulations in the development, con- 
struction, testing, and operation of X-ray optics has long been recognized. First 
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described by Wolter in 1952 [2!, the X-ray optics that would later be given his 
name were quickly adopted as the preferred configuration for X-ray imaging sys- 
tems [s]. In 1972, Chase and Van Speybroeck carried out the first numerical 
evaluation of the optical properties of Wolter optics. Since then several authors 
have investigated X-ray imaging systems (see for instance Q and references 
therein). 

To a large extent X-ray imaging systems can be treated as any other optical 
system, and their properties in terms of aberrations and focusing capabilities 
can be examined with the standard tools of optics. 

However, when detailed performance predictions of the effective area and the 
point spread function are required, specific ray-tracing tools are needed in order 
to take into account the change of the reflectivity with the angle of incidence, 
and the scattering caused by the details of the surface finish. To the best of 
our knowledge, and certainly within the domain of X-ray optics for astrophysical 
missions, this problem has always been solved with ad hoc tools built to address 
issues related to the particular design under consideration. 

We have developed a set of extensions for the Gcant4 toolkit that allows one 
to trace X-ray optics of arbitrary complexity, as long as the geometrical model 
of the system can be described by the geometry library available in Geant4. 
The core component of these extensions consists of a Geant4-compatible im- 
plementation of the reflection of X-ray photons on a surface. Additionally, by 
making use of the available Geant4 functionality, we have created a framework 
that allows one to introduce a microscopic description of the properties of the 
surfaces of the optics. Through this framework it becomes possible to model the 
effects that these microscopic properties have on the scattering of the photons. 

Besides introducing the notion of a generic tracer for X-ray optics, this de- 
velopment has the potential to ease the interaction of scientists and engineers 
during the study phase of X-ray space missions. Geant4-based models are in 
fact already widely used to study the effect of cosmic radiation on the spacecraft 
structures and instruments. By making use of the same tools in order to model 
both the spacecraft and the telescope, the study of a mission can be carried out 
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on the basis of a shared understanding of its mechanical and optical properties. 

This article is organized as follows. In §[5] we briefly describe the physics 
of X-ray reflection and scattering at grazing angles. In §[3] the X-ray oriented 
extensions to the Geant4 toolkit that we have implemented are described. We 
proceed then to describe a sample application of our extensions in §|4l Finally, 
in §[5] we draw some conclusions and make some final remarks. 

2. The Reflection of X-Rays 

The reflection of X-rays at the boundary between vacuum and a medium 
(see Fig. [1]) is described by the standard Fresnel equations (see for instance Q). 
Here, in keeping with the usage in the field, we use energy instead of wavelength. 

For a photon with grazing incidence angle Oi and energy E, the reflection 
coefficient R and the transmission coefficient T (both for the orthogonal and 
parallel polarizations) are: 



where n{E) is the medium's index of refraction. For different materials n{E) 
can be obtained from various sources. In practice, the reflectivity is calculated 
using experimentally determined optical constants that are valid in the X-ray 
regime. However, when these optical constants are not available because of the 
lack of measurements (in the case of non standard materials) or because the 
photon energy is in the gamma-ray regime, where only a few optical constants 
are known, one can construct the optical constants using the atomic form factors 
(e.g. B). 

The reflection coefficient given in Eq.(IT]) is valid in the case of an ideal, 
perfectly smooth surface. For a real surface, some fraction of the reflected 
photons will be scattered away from the specular direction. Under assumptions 
generally true for X-ray optics, it is customary to characterize a surface by its 
micro-roughness a, and decrease the reflectivity by the factor 
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In the context of astrophysical applications of X-ray optics, Eq.® is suffi- 
cient to perform calculations related to the effective area of a particular optical 
design. More complex calculations aiming at the prediction of the point spread 



function must use a more extensive theory (for instance [8|, whose details 
are affected by the type of surface data available (e.g., surface interferograms, 
profilometry measurements, atomic force telescope data). Our software im- 
plementation provides a generic interface for the application of any scattering 
model, as is described in the following section. 

3. X-Ray-Oriented Extensions of the Geant4 Toolkit 

The extension of the Geant4 toolkit with the physics of grazing angle scat- 
tering is realized through the implementation of three classes: 

G4XrayRefractionIndex This is an auxiliary class that manages the refrac- 
tion index data for a particular material. 
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G4XraySurface As mentioned in the previous section, the details of the in- 
teraction of X-rays on a surface can be described by scattering models 
of various complexity, and are also driven by the type of data describ- 
ing the microscopic surface details. It is therefore not surprising that 
even within the same model different surfaces may have different X-ray 
behaviours: for instance, an X-ray mirror will have one reflecting side, 
and a back side treated so as to absorb or scatter away from the main 
beam all photons that might fall on it. Also, when assembling an optical 
system for which disparate surface data types are available, it may be 
desirable to be able to assign different scattering models to different sur- 
faces. The class G4XraySurface provides a generic interface to do exactly 
this. It can be used with the standard Geant4 mechanisms to deflne ei- 
ther a logical boundary surface, or a logical skin surface. We have written 
sample classes that implement Presnel reflection as described in the previ- 
ous section, and that perturb the surface normal at the interaction point 
according to a Gaussian distribution. In general, G4XraySurface can be 
used to implement any scattering model appropriate for the particular 
situation considered. 

G4XrayGr£izingAngleScattering This class models a new Geant4 boundary 
process, and is applicable to any photon (the Geant4 particle gamma: 
no new particle needs to be defined). The process can compete with 
any other Geant4 processes that might be defined in the simulation (for 
instance photo-electric effect, or Compton scattering). When the Geant4 
tracking system detects that a photon is crossing the boundary between 
two volumes (normally from vacuum to a material), this process is invoked. 
In turn, the process queries the surface being crossed to see if it is an X- 
ray surface: if this is the case, the process delegates the actual interaction 
details to the surface. The scattered/reflected photon is handed back 
to the process who decides its fate — absorbed or transmitted, reflected, or 
propagated further — , and eventually hands the control back to the Geant4 
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tracking manager. The process can be activated so as to log the details 
of every X-ray interaction: the coordinates of the interaction point, the 
volume on which the interaction has taken place, the (possibly modified) 
local surface normal, incoming and outgoing directions of the photon, and 
the reflectivity used. This information, that is commonly used in the 
field to understand the details of the geometrical model and how certain 
features on the focal plane are formed, can be accessed using the standard 
track information facilities of Geant4. 

Using the extension a simple reflectometry set up can be constructed con- 
sisting of a single silicon plate as shown in Fig. HJa). The reflectivity for a 
fixed angle of incidence can then be determined, and compared with the data 
obtained using the tools made available by the Center for X-Ray Optics [l^, a 
standard reference in the field. In Fig. [DJb) the two data sets are shown to be 
in very good agreement. 



6 



(a) (b) 

Figure 3: (a) Photograph of the front view of the XMM-Newton telescope (courtesy of the 
European Space Agency), (b) The Geant4 model used for the comparisons. 



4. Application Example 

As an example we have modeled the X-ray optical system of the XMM- 
Newton telescope. Each of the three mirrors aboard XMM-Newton comprises 
58 concentric mirror shells in a Wolter-I configuration, where each mirror shell 
consists of polished, confocal parabolic and hyperbolic surfaces. The parabolic 
and hyperbolic solids that are included in the Geant4 package are not directly 
useful to construct a Wolter-I system, because they are filled solids. We therefore 
approximated them with many consecutive short conical shell segments with the 
appropriate slope. The front view of the telescope is shown in Fig.[3l^a). In this 
photograph the nest of mirror shells can clearly be seen together with the spider- 
like mechanical support structure. The total aperture ranges from 150 mm to 
350 mm. In Fig.[3ljb) the same geometry is shown as implemented in our model. 
The mirror shells are coated with a gold layer and we assume a surface roughness 
of 0.5 nm root-mean-square. 

In Fig. m we show the focal spot for an on-axis X-ray beam in comparison to 
the focal spot for an off-axis beam. The spot for the off-axis beam lies away from 
the center of the image plane as expected. For an off-axis angle of 2 • 10^^ degrees 
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Figure 4: (a) Focal spot of the XMM-Newton telescope for an on-axis point source, (b) Focal 
spot for an off-axis (0.02 degrees) point source. 



and a focal length of 7500 mm, the x-coordinate is expected to be 2.618 mm, in 
agreement with the simulated result of 2.622 mm. Note further that the shape 
of the spot becomes asymmetric due to coma, as expected. 

By Monte Carlo simulation we have determined the effective area of the 
XMM-Newton optics, by taking the blocking effect due to the mechanical sup- 
port structure into account. We have not taken into account any other effects, 
like dust contamination, or the exact form of the mirror shells. In Fig. [3 w e 
compare our results with the ground calibration measurements reported in 
our results are in very good agreement with the data, differing no more than 
5% from them. 



5. Conclusions 

We have created an extension of the Geant4 toolkit that makes it suitable 
as a platform with which to realize generic X-ray tracers. Besides modelling 
the grazing angle reflection of X-rays on surfaces, our extension can be used to 
model the effects of surface finish on the scattering of radiation. 

We have demonstrated the new capabilities of the toolkit by modelling one 
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energy [keV] 

Figure 5: On-axis effective area of the XMM-Newton telescope. The sohd Une shows the 
results of our simulation. The three data points (triangles) are taken from [ill . The difference 
between measured and simulated effective area is less than 5%. 
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of the Wolter-I mirrors used on the European Space Agency's X-ray observatory 
XMM-Newton. Without having access to the detailed metrological data that 
were used during the cahbration of XMM-Newton, we are able to reproduce 
the measured effective area to about 5%. We have also shown that the main 
properties of the point spread function are well reproduced. Being aware of the 
efforts that went into accurately modelling the response of the XMM-Newton 



mirrors 



12l | , we consider our result as a demonstration of the potential our work 
has to ease the design and calibration of future X-ray optics. 

With this extension we have turned Geant4 into a versatile modeler for X-ray 
optics. 
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Schematic representation of tlie reflection and transmission at the 

boundary between vacuum and material 

(a) Simulated reflectometry setup, (b) The reflectivity deter- 
mined from the simulation is compared with data retrieved from 
the Center for X-ray Optics [10 1 



(a) Photograph of the front view of the XMM-Newton telescope 
(courtesy of the European Space Agency) . (b) The Geant4 model 

used for the comparisons 

(a) Focal spot of the XMM-Newton telescope for an on-axis point 
source, (b) Focal spot for an off-axis (0.02 degrees) point source. 
On-axis effective area of the XMM-Newton telescope. The solid 
line shows the results of our simulation. The three data points 
(triangles) are taken from ll| . The difference between measured 
and simulated effective area is less than 5% 
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